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Abstract. We study a model of Brownian particles which are pumped with energy by means of a 
non-linear friction function, for which different types are discussed. A suitable expression for a 
non-linear, velocity-dependent friction function is derived by considering an internal energy depot 
of the Brownian particles. In this case, the friction function describes the pumping of energy in the 
range of small velocities, while in the range of large velocities the known limit of dissipative friction 
is reached. In order to investigate the influence of additional energy supply, we discuss the velocity 
distribution function for different cases. Analytical solutions of the corresponding Fokker-Planck 
equation in 2d are presented and compared with computer simulations. Different to the case of 
passive Brownian motion, we find several new features of the dynamics, such as the formation of 
limit cycles in the four- dimensional phase-space, a large mean squared displacement which increases 
quadratically with the energy supply, or non- equilibrium velocity distributions with crater-like form. 
Further, we point to some generalizations and possible applications of the model. 

PACS. 05.40.-a Fluctuation phenomena, random processes, noise, and Brownian motion - 05.45.- 
a Nonlinear dynamics and nonlinear dynamical systems - 82.20. Mj Nonequilibrium kinetics - 
87.15.Vv Diffusion - 87.15.Ya Fluctuations 



1 Introduction 

As we know from classical physics, Brownian motion de- 
notes the erratic motion of a small, but larger than molec- 
ular, particle in a surrounding medium, e.g. a gas or a liq- 
uid. This erratic motion results from the random impacts 
between the atoms or molecules of the medium and the 
(Brownian) particle, which cause changes in the direction 
and the amount of its velocity, v. 

This type of motion would be rather considered as pas- 
sive motion, simply because the Brownian particle does 
not play an active part in this motion. The Brownian 
particle keeps moving because the dissipation of energy 
caused by friction is compensated by the stochastic force, 
as expressed in the fluctuation-dissipation theorem (Ein- 
stein relation). 

In this paper, we are particularly interested in the 
non-equilibrium motion of Brownian particles. The out-of- 
equilibrium state is reached by considering an additional 
influx of energy which is transfered into kinetic energy 
of the particles. This will be considered by a more com- 
plex friction coefficient, 70, which now can be a space- 
and/or velocity-dependent function, j(r, v), which can be 
also negative. 
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Negative friction is known i.e. from technical construc- 
tions, where moving parts cause a loss of energy, which is 
from time to time compensated by the pumping of me- 
chanical energy. For example, in a mechanical clock the 
dissipation of energy by moving parts is compensated by 
a heavy weight. The weight stores potential energy, which 
is gradually transferred to the clock via friction, i.e. the 
strip with the weight is pulled down. Another example are 
violin strings, where friction-pumped oscillations occur, if 
the violin bow transfers energy to the string via friction. 
Already in the classical work "The Theory of Sound" of 
Lord Rayleigh considered motion with energy supply jl] . 

Provided a supercritical influx of energy, passive mo- 
tion could be transformed into active motion, which relies 
on the supply of energy from the surroundings. Active mo- 
tion is of interest for the dynamics of driven systems, such 
as physico-chemical ||-|(| or biological (7) systems. How- 
ever, recent models on self-driven particles often neglect 
the energetic aspects of active motion while focusing on 
the interaction of particles 

In order to describe active motion in the presence of 
stochastic forces, we have suggested a model of active 
Brownian particles, which have the ability to take up en- 
ergy from the environment, to store it in an internal energy 
depot and to convert internal into kinetic energy [jl0]-jl4j|. 
Two of us have shown that these particles could move in 
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an "high-velocity" or active mode, which results in differ- 
ent types of complex motion. 

Other previous versions of active Brownian particle 
models consider also specific activities, such as environ- 
mental changes or signal- response behavior. In particu- 
lar, the active Brownian particles are assumed to gener- 
ate a self-consistent field which effects their further move- 
ment 0-0. The non-linear feedback between the parti- 
cles and the field generated by themselves may result in 
an interactive structure formation process on the macro- 
scopic level. Hence, these models have been used to sim- 
ulate a broad variety of pattern formations in complex 
systems, ranging from physical to biological and social 
systems 

The main objective of this work is to study the in- 
fluence of different types of non-linear friction terms on 
active Brownian particles and the corresponding probabil- 
ity distributions. Hence, we neglect possible environmental 
changes of the particles and rather focus on the energetic 
aspects of motion in the presence of noise. In Sect.|, we in- 
troduce the idea of pumping by active friction and outline 
the basic dynamics of our model by means of Langevin and 
Fokker-Planck equations. Further, we investigate the dis- 
tribution function for spatially uniform systems, i.e. in the 
absence of an external potential. In Sect. || we consider ac- 
tive Brownian motion in external potentials. For the case 
of a parabolic potential, we derive analytical expressions 
for the distribution function and present computer simu- 
lations based on the Langevin equation, which show the 
occurrence of limit cycles. The basic ideas are then gen- 
eralized for other types of potentials. Finally, we discuss 
some possible applications for the dynamics of pumped 
Brownian particles. 



2 Active Brownian Motion of Free Particles 

2.1 Equations of Motion for Active Brownian Motion 

The motion of a Brownian particle with mass m, position 
r, and velocity v moving in a space-dependent potential, 
U{r), can be described by the following Langevin equa- 
tion: 



-Mv)v VU{r)+F(t) 

m 



(1) 



Here, j(v) is the friction coefficient, which is assumed to 
depend on velocity and and therefore implicitly on time. 
This rather general ansatz will allow us to consider more 
complex functions for the friction coefficient, as discussed 
below. J-(t) is a stochastic force with strength D and a 
5-correlated time dependence 



(jF(t)) = ; {T{t)T{t')) = 2D S(t - t') 



(2) 



In the case of thermal equilibrium systems with ^(v) = 
7o = const, we may assume that the loss of energy result- 
ing from friction, and the gain of energy resulting from 
the stochastic force, are compensated in the average. In 



this case the fluctuation-dissipation theorem (Einstein re- 
lation) applies: 

-£> = 7o (3) 

m 

T is the temperature and fee is the Boltzmann constant. 

In this paper we are mainly interested in the discussion 
of the probability density P(r, v, t) to find the particle at 
location r with velocity v at time t. As well known, the 
distribution function P(r,v,t) which corresponds to the 
Langevin Eq. ([!]), can be described by a Fokker-Planck 
equation of the form: 

8P(r,v,t) d ( , . „, . „8P(r,v,t) 

_ v dP(r,v,t) + L 8P(rv,t) 
or m ov 

In the special case -f(v) = 70 the stationary solution of Eq. 
(Q), P°(r,v), is known to be the Boltzmann distribution: 

p. (r ,„ ) = Cexp (_^M)^(_i^) (5) 

where the constant C results from the normalization con- 
dition. 

The major question discussed throughout this paper 
is how this known picture changes if we add a new de- 
gree of freedom to the model by considering that Brown- 
ian particles can be also pumped with energy from the 
environment. This extension leads to the model of active 
Brownian particles [|j^|ll],[T^,[l6| . While the usual dynamic 
approach to Brownian motion is based on the assumption 
of a passive friction described by a constant, non-negative 
friction coefficient, 70 > 0, active Brownian motion con- 
siders active friction as an additional mechanism of the 
particle to gain energy. In the present model, this will 
be described by a more complex friction coefficient "f(v) 
which is a velocity-dependent function. Several aspects of 
an additional space dependence were discussed in an ear- 
lier work ||. 

In the next section, we will introduce two possible 
ansatzes for the pumping of energy, where the active fric- 
tion depends on the velocity v of the particle. In Sect. || 
we will restrict the discussion to the case U(r) — const., 
but in Sect. || also the influence of external forces will be 
discussed. Further, from now on we will choose units in 
which m = 1. 



2.2 Pumping By Velocity-Dependent Friction 

The velocity-dependent pumping of energy plays an im- 
portant role e.g. in models of the theory of sound devel- 
oped by Rayleigh |jj already at the end of the last century. 
In the simplest case we may assume the following friction 
function for the Brownian particle: 



7(v) = -71 + 72 v = 71 — - 1 



(6) 
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This Rayleigh-type model is a standard model for self- 
sustained oscillations studied in different papers on non- 
linear dynamics [Bh,B3|]. We note that 



21 

72 



(7) 



defines a special value of the velocity where the friction 
function, Eq. (H), is zero (cf. Fig. m. 
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Fig. 1. Two different types of velocity dependent friction 
functions: (top) Rayleigh type, Eq. (^|), (bottom) Eq. (^). Pa- 
rameters: 71 = 5, 72 = 1.4, 70 = 2, vo = 1.889. 



Another standard model for active Brownian dynamics 
with a zero i>o in the friction function introduced in [Q 
reads (cf. Fig. |l|): 



70) = 7o (l - ^) 



(8) 



It has been shown that Eq. (|Sj) allows to describe the ac- 
tive motion of different cell types, such as granulocytes 
jj],[24|,^5| monocytes Q or neural crest cells ^7j. Here, 
the speed v expresses the fact that the motion of cells is 
not only driven by stochastic forces, instead cells are also 
capable of active motion. 

To compare Eqs. (JsJ) , (||) (cf. Fig. we note that in 
both ansatzes we have a certain range of (small) veloci- 
ties, where the friction coefficient can be negative, i.e. the 



motion of the particle can be pumped with energy. Due 
to the pumping mechanism, the conservation of energy 
clearly does not hold for the particle, i.e. we now have a 
strongly non-equilibrium canonical system ]2Sj] . 

Above a critical value of the velocity, described by Vq, 
the active friction turns into passive friction again. In Eq. 
(^]) , the increase of the friction with increasing velocity is 
not bound to a maximum value. However, this is the case 
in Eq. (|§|), where the friction function converges into the 
constant of passive friction, 70, for large v. The disadvan- 
tage of Eq. (g), on the other hand, is the singularity of the 
friction function for v — > 0. In Sect. 2.3, we will suggest 
a model of active friction which avoids these drawbacks, 
while providing a more general approach to derive a suit- 
able function 7(f). 

But before, we want to discuss some features which 
result from the Rayleigh ansatz, Eq. The stationary 
solution of the corresponding Fokker-Planck Eq. (|J) reads 
in the isotropic case: 



(9) 



The 



distribution P°(v), Eq. (|^), has its maximum for 
v = 0, if 71 < 0. For 71 > 0, the mean of the veloc- 
ity distribution is still zero, but the maxima are differ- 
ent. Considering a two-dimension motion, the distribution 
then has the form of a hat (cf. also Fig. ||), indicating that 
the particle most likely moves with a constant absolute 
velocity v 2 = Vq =71/72. 

For the calculation of the normalization constant C in 
Eq. (H) we use a method described by j?9) and find with 
72 = 1 the explicit expression pCj: 



c- 



2L 

AD 



1 + erf 



7i 
2-/D 



(10) 



Using Eqs. (g), ([l0|), we are able to calculate the different 
moments of the stationary velocity distribution and find: 



, 2M (2D) n 8 n „ . 

> - hk Wi Z2d(71 



with the function 

^2d(7l 



7tv ttU exp — — 
1 AD 



1 +erf 



7i 
2-/D 



(11) 



(12) 



In particular, the second moment, which is proportional 
to the temperature, and the fourth moment which is pro- 
portional to the fluctuations of the temperature, read: 



(v 2 ) = 7l +J- 

V 7T i 

1 + erf 

( W 4 )=7 1 2 + 2^ + 7i<^ 2 ) 



7i 
2VD 



(13) 



We note that a similar discussion can be carried out also 
for the case of the friction function, Eq. rf§). In this case 
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the stationary solution of the Fokker-Planck Eq. 
particular simplicity 0: 



P°(w) 



n f To / v 
XP \2D ( V ~ V °> 



is of 



(14) 



In the following section, we will derive a model which al- 
lows us to generalize both special cases introduced here, 
so that the further discussion can be unified. 



2.3 Pumping from an Internal Energy Depot 

In order to avoid the drawbacks of the friction functions, 
Eqs. (||), (||), mentioned in the previous section, we need a 
friction function 7(1?) which does not diverge in the limit of 
small or large velocities. Further the case of the "normal" 
passive friction should be included as a limiting case. The 
ansatz suggested here is based on the recently developed 



(i) 
(ii) 

(iii) 



model of Brownian motion with internal energy depot 1 10 

We assume that the Brownian particle should be able 
to take up external energy, which can be stored in an 
internal energy depot, e(t). This energy depot may be may 
be altered by three different processes: 

take-up of energy from the environment; where q is the 
pump rate of energy, 

internal dissipation, which is assumed to be propor- 
tional to the internal energy. Here the rate of energy 
loss, c, is assumed to be constant, 
conversion of internal energy into motion, where d (i> 2 ) 
is the rate of conversion of internal to kinetic degrees of 
freedom and r\ (i> 2 ) is the efficiency of conversion. This 
means that the depot energy may be used to accelerate 
motion of the particle. 

This extension of the model is motivated by investigations 
of active biological motion, which relies on the supply of 
energy. This energy can be dissipated by metabolic pro- 
cesses, but can be also converted into kinetic energy. The 
resulting balance equation for the internal energy depot 
of a pumped Brownian particle is then given by: 

j t e(t) = q - c e(t) - r,(v 2 )d(v 2 ) e(t) (15) 

Generalizing the ansatz given in Jn],[ll| we postulate the 
following Langevin equation for the motion of Brownian 
particles with internal energy depot: 

v + 7 f + VU = \ y ' e(t)v + T{t) (16) 

Let us consider now two limiting cases: 
(i) First, we assume 



q(r) = q ; d (v 2 ) = d 2 v 2 ; 77 (v 2 ) = 1 



(17) 



where d% > is the conversion factor. Using the assump- 
tion of a fast relaxation of the internal energy depot we 
find the quasistationary value: 



(18) 



With the quasistationary approximation eo, Eq. (fig), the 
internal energy can be adiabatically eliminated, and we 
find the following effective friction function: 



l{v 2 ) 



7o 



qo d 2 



(19) 



which is plotted in Fig. |. We see that in the range of 
small velocities pumping due to negative friction occurs, as 
an additional source of energy for the Brownian particle. 
Hence, slow particles are accelerated, while the motion of 
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Fig. 2. Velocity-dependent friction function 7(w), Eq. (^) vs. 
velocity v. The velocity ranges for "pumping" {^(v) < 0) and 
"dissipation" ("/(v) > 0) are indicated. Parameters: go — 10; 
c = 1.0; 70 = 2, d 2 = 0.7. 



The friction function, Eq. (|l9|), has again a zero, which 
reads in the considered case: 



<Zo 

7o 



c 

d~2 



(20) 



Eq. ( |T9| ) agrees with the Rayleigh ansatz, Eq. (H), in the 
limit of rather small velocities, v < Vq. From the condi- 
tions "f(v = 0) and j(v) — 0, we find for the two parame- 
ters of the Rayleigh ansatz: 



7i = 7o - 90- 



rf 2 



72 = 7o- 



d- 2 



(21) 



Thus, the further discussion will also cover this case as 
long as v < vq. 

(ii) In the second limiting case we consider a very large 
energy depot which is then assumed to be constant, e(t) = 
eo = const. However, we may assume that the efficiency 
of conversion, rj (i> 2 ) , decreases with increasing velocity as 
follows: 

.2\ _ Vi 



1] \ v 



1 



(22) 



where r\\ and 772 are constants. With 770 = ^1^2, this case 
leads to the effective friction function |13|: 



j(v 2 ) 



: 7o 



1 + 772V 2 



(23) 
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Both expressions, Eq. ( |19[ ) and Eq. (23), can be brought 
into the same form by defining some constants: 



7 (v 2 ) = 70 

Q = 

c7o 



1 - 



q 



1 + d*v 2 



7o 



; d* 



V2 



(24) 
(25) 



Using further the expression for the stationary velocity 
vl, Eq. ©, we find for Eq. © 



7 {v 2 ) = 70 



(<7o/7o) + (v 2 - v 2 ) 



(26) 



We restrict the further discussion to the two-dimensional 
space r = {x\, x 2 }. Then the stationary velocities v , Eqs. 
(Q), (pp[), where the friction is just compensated by the 
energy supply, define a cylinder in the four-dimensional 
space: 

v 2 +v 2 = v 2 (27) 

which attracts all deterministic trajectories of the dynamic 
system. The stationary solution of the Fokker-Planck Eq. 
(f|) reads for the friction function, Eq. (|l9|), and in the 
absence of an external potential, i.e. U(x\ 1 X2) = 0: 



p» = c(i + 



2\ 2D 



exp 



Compared to Eq. (|J) , which describes the Maxwellian ve- 
locity distribution of "simple" Brownian particles, a new 
prefactor appears now in Eq. ( p8| ) which results from the 
internal energy depot. In the range of small values of v 2 , 
the prefactor can be expressed by a power series truncated 
after the first order, and Eq. ( j28[ ) reads then: 



P°(v) 



exp 



7o 
2D 



1 - 



90^2 

C70 



(29) 



In Eq. (|29j) , the sign of the expression in the exponent de- 
pends significantly on the parameters which describe the 
balance of the energy depot. For a subcritical pumping of 
energy, qo^2 < C70, the expression in the exponent is nega- 
tive and an unimodal velocity distribution results, centered 
around the maximum vo = 0. This is the case of the "low 
velocity" or passive mode for the stationary motion which 
here corresponds to the Maxwellian velocity distribution. 
However, for supercritical pumping, 



q d 2 > C70 



(30) 



the exponent in Eq. (|29j) becomes positive, and a crater- 
like velocity distribution results. This is also shown in Fig. 
||. The maxima of P°(v) correspond to the solutions for 
Vq, Eq . ( 20j) . The corresponding "high velocity" or active 
mode [|12 [for the stationary motion is described by strong 
deviations from the Maxwell distribution. 
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Other parameters: 70 = 2, D = 2, c = 1, go = iu. mooe inai 
d 2 = 0.2 is the bifurcation point for the given set of parameters. 

2.4 Investigation of two Limit Cases 

In the limit of strong noise D ~ T — * 00, i.e. at high 
temperatures, we get from Eq. (pS) the known Maxwellian 
distribution by means of Eq. (pH: 



P°(v) = 



1 



2nk B T 



exp 



2k B T 



(31) 



which corresponds to standard Brownian motion in two 
dimensions. Hence, many characteristic quantities are ex- 
plicitly known as e.g. the dispersion of the velocities 

<V ) = 2k B T (32) 

and the Einstein relation for mean squared displacement 

k B T 



(r(t)-r(0) 



= 4- 



■i 



7n 



(33) 



In the other limiting case of strong activation, i.e. rela- 
tively weak noise D ~ T — > and/or strong pumping, we 
find a (^-distribution instead: 



P\v) = C5{v 2 -vl) ■ (v 2 )=v 2 



(34) 



In order to treat this case in the full phase space, we fol- 
low and introduce first an amplitude-phase represen- 
tation in the velocity-space: 



vi = v cos (if) ; v 2 = v sin (ip) 



(35) 
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This allows us to separate the variables and we get a dis- 
tribution function of the form: 

P(xi,x 2 ,Vi,V2,t) = P(x ll x 2 ,t) ■ 5(v\ + v 2 - Vq) ■ P(tp,t) 

(36) 

The distribution of the phase <p satisfies the Fokker-Planck 
equation: 

3 d 2 

-P(<p,t)=D lp7 ^P(<p,t) (37) 



dtp 2 



By means of the known solution of Eq. ( p7| ) we get for the 
mean square: 



D v t ; 



D 



(38) 



where D v is the angular diffusion constant. By means of 
this, the mean squared spatial displacement (r 2 (t)) of the 
particle can be calculated according to Q as: 



D 2 



( 2Dt\ 


- 1 









(39) 



For times t ^> Vq/D, we find from Eq. ( |39"| ) the following 
expression for the effective spatial diffusion constant: 



D eS = 



D 



2 
D 



(Jo 

To 



c 

d~2 



(40) 



where vq, Eq. (EG), considers the additional pumping of 
energy resultingfrom the friction function, 7(f). Due to 
this additional pumping, we obtain a high sensitivity with 
respect to noise expressed in the scaling with (1/D). 



3 Active Brownian Motion in External 
Potentials 

3.1 Motion in a Parabolic Potential 

In the following, we discuss the motion in a two- 
dimensional external potential, U(xi, x 2 ), which results in 
additional forces on the pumped Brownian particle. The 
case of a constant external force was experimentally and 
theoretically investigated in JtJ . Here we concentrate first 
on a simple non-linear potential: 



U(x 1 ,x 2 ) = -a{x 1 + x l) 



(41) 



If we for the moment restrict the discussion to a deter- 
ministic motion, the dynamics will be described by four 
coupled first-order differential equations: 



x\=Vi vi — -7 (vi, v 2 ) v\ - ax 1 
x 2 = v 2 v 2 — -7 (vi, v 2 ) v 2 - ax 2 



(42) 



For the one-dimensional Rayleigh model it is well known 
that this system processes a limit cycle corresponding to 
sustained oscillations with the energy En — 71/72- For the 
two-dimensional case we have shown |1 if that limit cycles 



occur in the four-dimensional space. The projection of this 
periodic motion to the {y\, v 2 } plane is again the cylinder, 
defined by Eq. (|27|). The projection to the {x\,x 2 } plane 
corresponds to a circle 

x\ + x\ = Tq = const. (43) 

The energy for motions on the limit cycle is: 

Eo = \ (vf + v\) + I (xf + 4) = \ 2 + %l (44) 

In p"l[ we have shown that any initial value of the energy 
converges (at least in the limit of strong pumping) to 



H 



En 



(45) 



This corresponds to an equal distribution between kinetic 
and potential energy. As for the harmonic oscillator in one 
dimension, both parts contribute the same amount to the 
total energy. This result was obtained in |ll| using the 
assumption that the energy is a slow (adiabatic) variable 
which allows a phase average with respect to the phases 
of the rotation. 

In explicite form we may represent the motion on the 
limit cycle in the four-dimensional space by the four equa- 
tions: 



xi = ?'o cos(ujt + tpa) 
x 2 = r sm(uit + ipo) 



i>i = — ?*o to sin(wt + ipo) 
v 2 = tq uj cos(ujt + (fo) 



(46) 

The frequency to is given by the time the particle needs 
for one period while moving on the circle with radius vq 
with constant speed vq. This leads to the relation 



too = — 
ro 



(47) 



This means that even in the case of strong pumping the 
particle oscillates with the frequency given by the linear 
oscillator frequency to. 

The trajectory defined by the four equations ( J46| ) is 
like a hoop in the four-dimensional space. Therefore, most 
projections to the two-dimensional subspaces are circles or 
ellipses. However there are two subspaces, namely {x\, v 2 } 
and {x 2 ,Vi}, where the projection is like a rod. 

A second limit cycle is obtained by time reversal: 



t 



-t; vi 



-v\] v 2 



-v 2 



(48) 



This limit cycle also forms a hula hoop which is different 
from the first one in that the projection to the {£1,2:2} 
plane has the opposite rotation direction compared to the 
first one. However both limit cycles have the same pro- 
jections to the {xi,x 2 } and to the {^1,^2} plane. The 
separatrix between the two attractor regions is given by 
the following plane in the four-dimensional space: 



Vi + v 2 = 



(49) 



Applying similar arguments to the stochastic problem we 
find that the two hoop-rings are converted into a dis- 
tribution looking like two embracing hoops with finite 
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size, which for strong noise converts into two embracing 
tires in the four-dimensional space E3. The projections 
of the distribution to the {xi, X2, ^2} subspace and to the 
{t>i, «2, £2} subspace are two-dimensional rings as shown 
in Fig. ||. While in the deterministic case either left- or 
righthanded rotations are found, in the stochastic case 
the system may switch randomly between the left- and 
righthand rotations, since the separatrix becomes pene- 
trating (see Fig. ||). This can be also confirmed by looking 
at the projections of the distribution to the {x±, V2} plane 
and the {X2 , v\ } plane (cf . Fig. ||) . Since the hula hoop 
distribution intersects perpendicular the {x\,V2\ and the 
{X2, v\} plane, the projections to these planes are rod- like 
and the intersection manifold with these planes consists 
of two ellipses located in the diagonals of the planes. 
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In order to get the explicite form of the distribution 
we may introduce the amplitude-phase representation 



x\ = g cos(ujt + ip) v\ — — quj sin(w£ + (p) 
X2 = Q sin(W + ip) t>2 = quj cos(u;i + p) 



(50) 



where the radius g is now a slow stochastic variable and 
the phase p is a fast stochastic variable. By using the 
standard procedure of averaging with respect to the fast 
phases we get for the Rayleigh model of pumping, Eq. (||), 
the following distribution of the radii: 



Fig. 4. Projections of a stochastic trajectory on the three- 
dimensional subspaces {xi, X2, V2} and {xi,X2, V2} for very low 
noise (D = 0.01). The particle will not leave the attractor 
which was reached dependant on the initial conditions (One 
of the trajectories represents righthanded and the other on 
lefthanded rotations). The trajectories were obtained with long 
time computer simulations of Eq. ( |l6| ) with constant energy 
depot, Eq. (|l8|). Other parameters: 70 = 0.2, c = 0.01, d,2 = 
0.1, a = 2.0 and q = 10.0. 



Po(q) - exp 



1-4 



(51) 
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This distribution has its maximum value at g = rrj. 

For the friction function Eq. (19) resulting from the 
depot model, we get the following stationary solution for 
the distribution of the radii: 



Po(q) 



«2 9 

1 + — ag 2 

c 



•10 

2 ; > 



exp 



70 2 N 

ag 

D y j 



(52) 
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Fig. 6. Poincare plots and projections of a stochastic trajec- 
tory on the two-dimensional subspaces {xi,V2} and {x2,vi} 
(same simulation as in Fig. |^). The noise strength is large 
enough to allow one particle to penetrate the separatrix be- 
tween the two hula-hoop-like attractors in finite time. For 
smaller noise one would need to simulate the particle's mo- 
tion with different initial conditions to reach both hoops/tires. 



Again, we find that the maximum is located at: 



Q = r o = — = 

0J Z 



(Jo 

To 



c 

T 2 



i 



(53) 



Due to the special form of the attractor which consists, 
as pointed out above, of two hula hoops embracing each 
other, the distribution in the phase space cannot be con- 
stant at g = const. An analytical expression for the dis- 
tribution in the four-dimensional phase space is not yet 
available. But it is expected that the probability density 
is concentrated around the two deterministic limit cycles. 



arg, is valid, which leads to the relation: loq — vo/ro = w. 
For other radially symmetric potentials U(r), this relation 
has to be replaced by the condition that on the limit cy- 
cle the attracting radial forces are in equilibrium with the 
centrifugal forces. This condition leads to 



? = \U'(r )\ 



(54) 



For a given Vq, this equation defines an implicit relation 
for the the equilibrium radius ro, namely: Vq = rg \ U'(ro)\. 
Then the frequency of the limit cycle oscillations is given 

*S-$-PM- (55) 

In the case of linear oscillators this leads to ojq — y/a as 
before, Eq. (f47j). But e.g. for the case of a quartic potential 



U(r) = ->, 

we get the limit cycle frequency 



LU 



fc 1 / 4 
1/2 



(56) 



(57) 



If Eq. ( |54| ) has several solutions for the equilibrium radius 
ro, the dynamics might be much more complicated, e.g. 
we could find Kepler-like orbits oscillating between the 
solutions for ro- In other words, we then find - in addition 
to the driven rotations already mentioned - also driven 
oscillations between the multiple solutions of Eq. ( |54| ) . 

In the general case of potentials U(xi,x 2 ) which do 
not obey a radial symmetry, the local curvature of po- 
tential levels U {x\ : x 2 ) = const, replaces the role of the 
radius ro- We thus define qjj (xi,x 2 ) as the radius of the 
local curvature. In general, the global dynamics will be 
very complicated, however the local dynamics can be de- 
scribed as follows: Most particles will move with constant 
kinetic energy 



v 2 = v 2 +v 2 =v 2 = const. 



(58) 



along the equipotential lines, U (x±, x 2 ) = const. In depen- 
dence on the initial conditions two nearly equally proba- 
ble directions of the trajectories, a counter-clockwise and 
a clockwise direction, are possible. Among the different 
possible trajectories the most stable will be the one which 
fulfills the condition that the potential forces and centrifu- 
gal forces are in equilibrium: 



Vq = QU (xi,X 2 ) |Vt/ (xi,X 2 )\ 



(59) 



This equation is a generalization of Eq. (B4|) . 



4 Discussion and Conclusion 



3.2 Generalization for Other Potentials 

For the harmonic potential discussed above the equal dis- 
tribution between potential and kinetic energy, mi>g = 



An interesting application of the theoretical results given 
above, is the following: Let us imagine a system of ac- 
tive Brownian particles which are pairwise bound by a po- 
tential U(ri — r 2 ) to dumb-bell-like configurations. Thus, 
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the pairs of active particles could be regarded as "active 
Brownian molecules" . Then the motion of each "molecule" 
consists of two independent parts: the free motion of the 
center of mass and the relative motion under the influence 
of the potential. 

The motion of the center of mass is described by the 
equations given in Sect. 2, while the relative motion is 
described by the equations given in Sect. 3. As a conse- 
quence, the center of mass of the dumb-bell will perform 
a driven Brownian motion, but in addition the dumb-bell 
is driven to rotate around the center of mass. We then 
observe that a system of pumped Brownian molecules - 
with respect to their center of mass velocities - has a dis- 
tribution corresponding to Eq. (S) or Eq. (|28|). However, 
since the internal degrees of freedom are excited, we also 
observe driven rotations and in general (if Eq. ( p4| ) has 
at least two solutions) also driven oscillations. Thus, we 
conclude that the mechanisms described here may be used 
also to excite the internal degrees of freedom of Brownian 
molecules. 

Another possible application is the motion of clusters 
of active Brownian molecules. Similar to the case of the 
dumb-bells, these clusters are driven by the take-up of en- 
ergy to perform spontaneous rotations. Eventually, a sta- 
tionary state will be reached which is a mixture of rotating 
clusters or droplets, respectively. 

These suggested applications for the dynamics of 
pumped Brownian particles will of course need further in- 
vestigations. In this paper the aim was rather to study the 
influence of non-linear friction terms on Brownian parti- 
cles from a more general perspective. 

A suitable expression for a non-linear, velocity- 
dependent friction function j(v) has been obtained by 
considering an internal energy depot of the Brownian par- 
ticles. This energy depot can be changed by energy take- 
up, internal dissipation and conversion of internal into ki- 
netic energy. Provided a fast relaxation, the depot can 
be described by a quasistationary value. We then found 
that the velocity-dependent friction function describes the 
pumping of energy in the range of small velocities, while 
in the range of large velocities the known limit of (dissi- 
pative) friction is reached. 

In this paper, the focus was on the distribution func- 
tion of pumped Brownian particles. Therefore, in addition 
to the mechanisms of pumping, we had to consider the de- 
pendence on stochastic influences, namely on the strength 
of the stochastic force, D. We could find new dynamical 
features for pumped Brownian particles, such as: 

(i) new diffusive properties with large mean squared dis- 
placement, which increases quadratically with the en- 
ergy supply and scales with the noise intensity as 1 / D 

(ii) non-equilibrium velocity distributions with crater-like 
shape, 

(iii) the formation of limit cycles corresponding to 
left/righthand rotations, 

(iv) the switch between these opposite directions, i.e. the 
penetration of the separation, during the stochastic 
motion of a pumped Brownian particle. 



Qualitative similar features are obtained for the motion of 
driven particles in physico-chemical systems, for example 
for the motion of liquid droplets on hot surfaces, or the 
motion of surface-active solid particles jjj-Jf]. There are 
also relations to active biological motion (7,[ll|]. In this 
paper, we did not intend to model any particular object 
but analyzed the general physical non-equilibrium proper- 
ties of such systems. While our investigations are based on 
rather simple physical assumptions on non-linear friction 
for Brownian particles, we found a rich dynamics, which 
might be of interest for more applied investigations later. 
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